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Number-resolving photo-detection is necessary for many quantum optics experiments, especially 
in the application of entangled state preparation. Several schemes have been proposed for approx- 
imating number-resolving photo-detection using non-number-resolving detectors. Such techniques 
include multi-port detection and time-division multiplexing. We provide a detailed analysis and 
comparison of different number-resolving detection schemes, with a view to creating a useful ref- 
erence for experimentalists. We show that the ideal architecture for projective measurements is 
a function of the detector's dark count and efficiency parameters. We also describe a process for 
selecting an appropriate topology given actual experimental component parameters. 



I. INTRODUCTION 

Number-resolving photo-detection is a necessary pre- 
requisite for many important quantum optics appli- 
cations, most notably in the emerging field of opti- 
cal quantum information processing [l|, Unfortu- 
nately, presently available photo-detectors are incapable 
of resolving photon number with high fidelity. In fact, 
most commonly available photo-detectors are so-called 
'bucket' or 'on/off' detectors, which can distinguish only 
between two cases - no photons, and one or more pho- 
tons. These practical limitations in the number-resolving 
capabilities of photo-detectors have motivated the de- 
velopment of techniques for approximating number- 
resolving detection using non-number-resolving detec- 
tors. Most notably, such techniques include multi-port 
networks [^, 3, IE Q , time-division multiplexing (TDM) 
0) B B [U EI and visible light photon counting mod- 
ules (VLPC's) i, [H, [11. All of these techniques are 
variations on a single simple idea - the optical field is 
distributed across multiple modes which are measured 
independently. For a large number of modes the proba- 
bility of any given mode being populated by more than 
one photon approaches zero. Thus, the sum of the de- 
tection events across the modes closely approximates the 
number of photons in the incident state. 

In this paper we analyze and compare different archi- 
tectures for implementing photon-number-resolving de- 
tection. Our analysis is primarily concerned with the 
performance of photon-number-resolving projective mea- 
surements. Such measurements are central to the prepa- 
ration of, for example, one mode of a bipartite entangled 
state by the detection of the other as described in dj|- 
Whilst photon number resolving detectors are also useful 



for the reconstruction of photostatistics, this application 
is not considered in this paper. Statistical reconstruction 
including the effects of the most dominant practical ef- 
fects - loss and dark-counts have been extensively studied 
in [3 and 

We begin by considering a very general setting, a lossy 
multi-port detector cascade i^. Most architectures for 
implementing number-resolving detection we are aware 
of reduce to specific instances of this setting. First we de- 
rive an expression for the conditional probability P{m\n) 
of detecting m photons given an incident 7i-photon Fock 
state. Our analysis includes the experimental effects of 
loss, finite detection eflaciency and dark counts neglected 
in other similar treatments [8| . We omit consideration of 
after-pulsing, as this may be disregarded with appropri- 
ate detector timing. We then derive a positive operator 
value measure (POVM) description of the measurement 
process, which is directly related to these conditional 
probabilities. The POVM description can be related to a 
quantum process description, which provides a very gen- 
eral characterization that can be applied to understand 
the dynamics of many different systems. 

Using this approach we characterize and contrast the 
A'^-port detector, and two primary variants of time divi- 
sion multiplexed (TDM) detectors - the balanced TDM 
detector JXl] and the loop detector from a projec- 
tive measurement perspective. Finally, we conclude with 
a design procedure for selecting an appropriate architec- 
ture and then optimising the design. We demonstrate 
that given contemporary experimental component pa- 
rameters, the balanced TDM architecture is generally 
most useful. 



II. GENERAL ANALYSIS 



We begin by considering a very general setting as fol- 
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photon Fock state incident upon one of the inputs (see 
Fig. [T]). The input couples to the outputs with couphng 
efficiency Pc{i) between the input and output i. That is, 
a given photon has probabihty pdi) of being directed to 
output i. The paths through the interferometer are lossy. 
The probability of a photon following the path to output 
i being lost is pioss(*)- Following the A^-port, each output 
is measured independently using a bucket detector [l9[. 
Each detector is subject to a dark-count probability of 
Pdc- We do not need to explicitly introduce a detector 
efficiency term, since this can be absorbed into the pioss 
terms. 
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FIG. 1: A''-port detector cascade. A state is incident upon 
one of the inputs, while vacuum states are input into the 
remaining inputs. The A'^-port then applies some transforma- 
tion which distributes the incident field across the outputs 
1 . . . N. Finally, the outputs are detected independently using 
non-number-resolving detectors. 



A. Measurement probability analysis 

Upon detection of the outputs of the cascade we will 
observe a particular pattern of detector clicks. We re- 
fer to a particular combination of clicks as a signature. 
Thus, a signature is simply a binary sequence of length 
N. Let us consider the probability of a particular signa- 
ture occurring. We use the notation Psig(rf) to denote the 
probability of the detection signature occurring whereby 
each of the detectors in vector d trigger, and all others 
do not. This probability is given by I20l| 
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The first line of this expression is a sum over all possible 
configurations of how the n photons can reach the N 
outputs (before loss) , weighted by the probability of each 
event. The variables rii are the number of photons that 
reach the i^^ output, before loss. The second and third 
lines represent the probabilities of the desired signature 
occurring, given the respective configuration. 

The net probability of detecting m photons is given by 
summing over all combinations of detection signatures 



where \d\ = m, where |d| is total number of clicks. Thus, 
P{m\n) = (2) 

\d\—m 

Note that a complete set of probabilities P{m\n)\/ m,n 
completely characterizes the operation of the detector. 

It is worth noting that Psig(rf) is not necessarily the 
same for all d where \d\ — m. That is, different signatures 
corresponding to the same measured number of photons 
needn't have equal probabilities of occurring. This will 
be the case, for example, in loop-based time-division- 
multiplexers, where there is an inherent asymmetry in 
photon arrival probabilities across different time-bins. 
This is something we will discuss in detail in later sec- 
tions. 

Qualitatively one can make several observations about 
any scheme which is an instance of this general scenario. 
First, in the ideal case (pioss — Pdc ~ 0) we expect the TV- 
port to approach an ideal detector in the limit N ^ oo 
and Pc{i) V z. In other words, in the limit where 
the incident field is split into an infinite number of in- 
finitely small components the probability of multiple pho- 
tons reaching a single detector approaches zero and the 
probability of correctly measuring photon number ap- 
proaches unity. For pdc > however this will no longer 
be the case, since each output port introduces an extra 
opportunity for a dark-count to occur. Thus, for large N 
false counts become a certainty. Therefore, for any given 
application one expects there to be an optimal value of 
N subsequently referred to as A'opt , which depends upon 
the experimental parameters. Based on this observation 
one intuitively expects that for small pdc it will be better 
to use large N, whereas for large pdc it will be better to 
use smaller N. 



B. POVM & quantum process descriptions 

We have derived a general expression for the opera- 
tion of a lossy A^-port detector in terms of the condi- 
tional probability P(m\n). For many applications a sim- 
ple expression for the conditional probability is insuffi- 
cient. Most notably, when using a detector to implement 
projective photon number measurements, a simple prob- 
ability measure is insufficient to derive the form of the 
projected state. For this a description of the measure- 
ment process in terms of measurement operators is nec- 
essary. We now consider the POVM description of the 
measurement process implemented by this general detec- 
tor. We denote the POVM element corresponding to the 
TO-photon detection outcome by II(m). This POVM ele- 
ment takes the form 



n(m) = YP{m\n)\n){n\ 



(3) 



Also, we implicitly assume that the measured state is 
traced out after measurement. This is because the mea- 
surement process destroys the incident state. Note that 
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for an ideal detector (i.e. N cxd, pc 0, pioss = Pdc = 
0), we have P{m\n) — Sm,n, in which case the POVM 
elements reduce to 

Ilidoail'Ti) |m)(m|, (4) 

as expected for an ideal number-resolving detector. In 
the general case however, these POVM's result in mix- 
ing over different possible measurement outcomes. This 
occurs for three independent reasons: loss means that an 
n-photon state may be confused for a < n-photon state; 
dark-counts mean that an n-photon state may be con- 
fused for a > n-photon state; and finite N means that 
the probability of more than one photon appearing at a 
given output is non-zero, again meaning that an n-photon 
state may be confused for a < n-photon state. 

The measurement process may also be described as a 
quantum process. When expressed in the basis of photon 
number projectors these processes are characterized by 
diagonal process matrices, 

= Y,PMn)EnpEr^ (5) 

n 

where x is a process matrix and £"„ = |n)(n| are 
the photon number projectors. In the ideal case where 
P{m\n) = 5m,m the process matrix corresponding to this 
quantum process will be the zero matrix with a single 
'1' at the to'^ location along the main diagonal. Thus, 
[x{'m)]i.j — Si.jSj,m and Eq. [5] reduces to 

£mip) = EmpErn (6) 

III. EXAMPLE DETECTION ARCHITECTURES 

In this section we apply our general analysis to sev- 
eral specific well-known architectures that have been 
experimentally demonstrated. We first consider the 
case of a balanced iV-port detector, perhaps the best 
known scheme for implementing number-resolving photo- 
detection. Then we consider two variations on time- 
division-multiplexed photo-detection. These two varia- 
tions are distinct in that their loss and coupling charac- 
teristics are inherently different. In one case the coupling 
terms are uniform, while in the other they are necessar- 
ily non-uniform. The relevant parameters for the different 
detection architectures considered are thus summarized 
in Table U 

A. Balanced A'^-port detection 

We first consider the case of a balanced iV-port detec- 
tor of Fig.[TJ Here an incident state is distributed equally 
across N outputs and thus the coupling parameters are 



Architecture pcji) pioss(i) 

Balanced A''-port 1/N pioss 

Loop TDM pcil-pcY-^ l-tltlt}-^r]d,t 
Balanced TDM 1/N 1 - t]rh]^+^r]det 

TABLE L Summary of relevant parameters for different de- 
tection architectures. The parameters ti, tc and ts indicate 
the transmission (i.e. 1— loss) of the optical fibre, coupler and 
switch components illustrated in Figs. [2] and [3] respectively, 
m is the number of stages of the balanced TDM detector. 
?7dct is the quantum efficiency of the final non photon number 
resolving detector used. 

all equal, Pc{i) — l/N V i. If it is assumed that the device 
is internally constructed such that loss is uniform across 
the outputs (i.e. by using a tree network of ideal beam- 
splitters), pioss(«) = Pioss V i. Such a symmetric device 
has been previously analysed in detail in 

In practise, due to their complexity, experimentalists 
do not construct TV-ports directly. Instead they build se- 
tups which closely approximate A^-port detection. Per- 
haps the closest approximation is the visible light photon 
counter (VLPC). Here an incident light field is spread 
out spatially and incident upon are large detector con- 
sisting of many small active detection areas which may 
trigger independently. Although this provides a close ap- 
proximation of a balanced A'^-port, it is imperfect due to 
asymmetry in the distribution of the field across the de- 
tector regions. Specifically, the light field incident upon 
the detector will have a roughly 2D Gaussian distribu- 
tion, rather than a uniform distribution. Nonetheless, in 
the regime where the number of detector regions is large 
compared to the number of incident photons this pro- 
vides a good approximation of number-resolving photo- 
detection. 



B. Loop time-division multiplexed detection 

Now we consider the loop-detector illustrated in Fig.[21 
This architecture has previously been considered in Ref. 
[ll[. Here the optical state is coupled into a fibre loop. 
The loop couples out again via a coupler with coupling 
strength pc- Thus, after each round-trip of the loop pho- 
tons have a probability pc of coupling out to the photo- 
detector. 

We relate the loop-detector to the generalized A^-port 
interferometer by letting the i*"^ port of the interferom- 
eter represent the corresponding time-bin of the loop- 
detector. Since photons in the i^^ (where i > 1) time-bin 
must undergo i — 1 round-trips and make i transits of the 
switch and coupler we havepioss(*) = 1 — is^c^}~^??dct and 

Pc{i) ^Pc{l~PcY~^- 

In principle a loop detector can continue detecting for 
an arbitrarily long time span, thereby effectively imple- 
menting N oo. In practise however, because of the ex- 
ponential decay in the coupling term with the number of 
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photo-detector 

FIG. 2: Loop time-division multiplexing. A state enters the 
fiber on the left and is coupled into the fiber loop by the 
electro-optic switch / coupler. It then repeatedly circulates 
through the fiber loop. Following each round-trip photons 
have some probability of coupling out to the photo-detector 
through the coupler. Thus, the state is divided into discrete 
time-bins, which are independently detected. t{ refers to the 
optical transmission of the fibre loop, tc is the transmission of 
the coupler, including connector losses, ts indicates the trans- 
mission through the switch. 

round-trips, it is sufRcient to truncate measurement and 
consider a relatively small number of time-bins. In fact, 
in the presence of dark-count it is preferable to truncate 
the number of measured time-bins. The reason for this 
is that the contribution of legitimate counts drops ex- 
ponentially against the number of round-trips, whereas 
the dark-count rate stays constant. Thus, for large N the 
later time-bins will achieve nothing other than to con- 
tribute unwanted dark-counts. 

Truncation gives rise to errors. Specifically, there will 
be some probability that photons land in the truncated 
region, which are discarded by the analysis, giving rise to 
erroneous results. In an experimental context this makes 
the measurement results ambiguous, reducing the fidelity 
of the measured state, while in a theoretical context this 
introduces an error margin in the analysis. Thus, one 
must be careful when considering results for the loop 
TDM to ensure that we are in a regime where the proba- 
bility of photons landing in the truncated region is small. 
Classically, this probability is of the order 

terror = (l-Pc)"^. (7) 

For quantum optical states this probability may be much 
lower, depending upon the likelihood of photon numbers 
> TV. 



C. Balanced time-division multiplexed detection 

The final architecture we consider is that of balanced 
time-division multiplexing, shown in Fig. [31 This archi- 
tecture consists of m stages (not to be confused with to, 
the number of measured photons - the context will make 
the distinction clear) , giving rise to = 2™ distinct non- 
overlapping time bins and is topologically equivalent to a 
balanced iV-port. That is, the coupling to each time bin 



is equal and Pc{i) = l/N V i. Despite being topologically 
equivalent, the balanced TDM differs from balanced N- 
port detection in that loss rates are not uniform across 
each time bin. Specifically, photons reaching the i*'' bin 
will pass through i — 1 lengths of fiber and to + I 50:50 
couplers. Hence we have pioss(i) = 1 — ty^t^c^^Vdct- 



2t, 




FIG. 3: A balanced TDM setup for m = 2 stages (i.e. iV = 
4). Thus, photons are distributed equally across four possible 
paths. t{ and are as defined previously, in Fig. [51 

While the iV-port and loop TDM architectures exhibit 
an optimal number of bins A'opt, which optimizes fidelity, 
and a minimum number of bins required to resolve the de- 
sired photon number, Nj^in, the balanced TDM architec- 
ture exhibits an optimal and minimum number of stages 
TOopt and TOinin respectively. 

IV. ANALYSIS OF ARCHITECTURES FOR 
THE CONDITIONAL PREPARATION OF FOCK 
STATES VIA PARAMETRIC 
DOWN-CONVERSION 

To introduce our experimental example and to expand 
on the general outline provided in the previous section, 
we now consider a common scenario; conditional prepara- 
tion of specific Fock states via non-degenerate paramet- 
ric down-conversion. This example is very important to 
present-day experiments where this is the defacto stan- 
dard for the preparation of single photons. Here we con- 
dition on detecting some number of photons in one arm 
of the down-converter. When conditioning succeeds we 
expect an equal number of photons to be present in the 
other output due to the photon number correlations be- 
tween the two output modes. 

The output state of a non-degenerate down-converter 
takes the form 

\^)^^=L=J2x''\n)\n), (8) 

where x is a parameter related to the down-conversion 
strength (not to be confused with the process matrix 
we introduced previously). Next we apply the to pho- 
ton measurement process, £m, to obtain the conditioned 

state, 

Pcond = £m{W{i>\) 

= -^^P(™|n)x^"|n)(n|. (9) 
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Following renormalization, the prepared state is given by 



N=2,Pdc=0 



Pprcp 



Pcond 



'cond ) 



X;„-P(m-|n)x^"|n)(r 
E„P(m|n)x2" 



(10) 



To quantify how well the state preparation procedure 
works, we calculate the fidelity between the the condi- 
tionally prepared state and the expected state \m), 



F{pcond, \m)) = 



2m 



P{'m\m)x 
X;„P(m|n)x2 



(11) 



In the limit of ideal photo-detection, P{m\n) = 5m,m the 
fidelity reduces to unity, as expected. Notice that that 
non-zero P{m\m), as is the case for non-ideal photo- 
detection, adversely affects the detection fidelity. The 
above fidelity only provide a measure of how sure we 
are of the desired state upon a detection event. The ac- 
tual probability of the desired (ml projection event pdet 
is given by 



Pdct 



= 5]P(m|n)|(n|^}|2 = ^^^P(mH 



(12) 

In the case of non-ideal detection, P(m\m) will in general 
be < 1 and the terms P{m\n) (m ^ n) will be non-zero. 



A. Balanced A'^-port detection 

First let us consider the operation of this setup using 
a balanced A^-port configuration for the conditioning de- 
tector. 

In Fig.|4]we plot the fidelity of projection onto a single 
photon state against x for different values of pioss, where 
we only assume detector loss, and pdc- As expected, the 
fidelity decreases monotonically with both these param- 
eters. In the absence of dark-counts the fidelity drops 
monotonically with x- This is because as x increases so 
too does the probability of higher photon number terms 
arising. In the presence of dark-counts, we observe max- 
ima in F against x- The initial increase in F is because 
dark-count rates are constant whereas the probability of 
a single photon number term occurring is monotonically 
increasing. Thus, when the probability of a single photon 
term occurring is very small the effects of dark-counts 
will dominate, reducing fidelity. As x increases so too 
does the effect of the desired single photon terms. Then 
for even higher values of x the probability of higher pho- 
ton number terms becomes significant, also reducing the 
fidelity. 

In Fig. [S] we plot the fidelity of the conditionally 
prepared state against down-conversion strength, x, for 
N — 2 and N ~ 5, and various values of pdc- In the 
ideal case, pdc — Pioss = 0, it is beneficial to use larger N 
since this increases the confidence when discriminating 
between different number states. However, in the pres- 
ence of dark-counts it is not strictly beneficial to use 
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Ploss =0.1 
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FIG. 4: Fidelity against loss (top) and dark-count (bottom) 
rates for fixed N = 2. 



larger N, since increasing N also increases the proba- 
bility of a dark-count occurring. 

This observation motivates us to determine iVopt for 
given loss and dark-count parameters. In Fig. [6] we plot 
the value of iVopt (i.e. the value of TV that optimized F) 
against x ^nd pdc- There are two important trends taking 
place. First, as the down-conversion strength increases, it 
becomes increasingly desirable to increase N . This is be- 
cause as X increases so too does the probability of gener- 
ating higher photon numbers, thus making the additional 
number resolving power of higher N necessary. Second, as 
the dark-count rate increases it becomes desirable to use 
smaller N. This is because the additional ports increase 
the probability of a dark-count occurring. We have not 
included a plot against pioss in this case, since it's effect 
on iVopt is close to uniform and almost negligible for the 
range of x £^nd pdc considered. 



B. Loop time-division multiplexing 

Next we consider the performance of loop TDM de- 
tection in the context of conditional photon preparation. 
The inherent asymmetry in the loop TDM architecture 
suggests there may also be asymmetry in the fidelity of 
different detection signatures corresponding to the same 
number of photons. For example, consider the case where 
there is no loss, but some dark-counts occurring. The 
dark-count probability is uniform across the detectors, 
but the probability distribution of incident photons is 
biased towards the earlier time-bins. Thus, a detection 
event occurring in an earlier time-bin will have higher rel- 
ative probability of being caused by an incident photon 
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Pdc = 0.001, Pioss = 
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FIG. 5: Comparison of the fidelity of a conditionally prepared 
state via non-degenerate parametric down-conversion for the 
balanced A^-port detection architecture. 



than if the detection had occurred at the second time-bin. 

To illustrate this, in Fig. [7] we plot the fidelity against 
Pdc and ti (and assume tc,ts = 1) for the {1, 0, 0, 0, 0} and 
{0, 1, 0, 0, 0} signatures, both of which correspond to de- 
tection of a single photon. In the presence of dark-counts 
we indeed observe that it is favorable to condition on the 
{1,0,...} signature than {0, 1, . . . }. In the presence of 
loss the converse is true. However, this does not arise be- 
cause of the loss itself. Rather it is inherent - higher order 
time bins have exponentially lower coupling probabilities. 




FIG. 6: Optimal value of A*' against down-conversion strength 
X, and dark-count rate pdc for a balanced A''-port configura- 
tion. Ploss = 0. 




0.01 0.02 0.03 0.04 0.05 

Pdc 




0.2 0.4 0.6 0.8 1 

Ploss 



FIG. 7: Asymmetry in the performance of different detection 
events corresponding to |d| = 1 in the loop TDM architecture. 
X = 0.1, pc = 0.5. 

and therefore stronger number resolving power. 

Thus, loss and dark-counts are competing parameters 
when determining the optimal conditioning signature. 
However, note that the dependence of fidelity on dark- 
counts is much greater than on loss. There are two rea- 
sons for this. First, there is a multiplicative effect from 
dark-counts - i.e. every detector is independently sub- 
ject to dark-counts. Thus, in our simulation where N = 5, 
there will be a roughly 5 times multiplicative factor in the 
dependence on dark-counts. This highlights the necessity 
of truncating loop TDM to minimize dark-count multi- 
plication. Second, loss causes n > 1 terms to be confused 
as n = 1 terms. However, because n > 1 terms have very 
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low probability, the fidelity does not suffer very much 
from these confusions. This is in contrast to the dark- 
count effect, which occur independently of how probable 
the n > 1 terms are. 



C. Balanced time-division-multiplexing 

Finally we consider operation when using balanced 
TDM detection. Qualitatively the behavior is very similar 
to the iV-port considered previously. In particular, in the 
presence of dark-counts there are distinct regions where 
different values of N are optimal. In the presence of no 
fiber loss (but potentially detector loss), we do not ob- 
serve that the fidelity varies significantly across different 
signatures corresponding to the same number of detected 
photons (the graphics for this have not been included). 
This arises because this system is almost balanced, so ide- 
ally the confidence is distributed evenly across the bins. 
There is a slight deviation from this because, although 
coupling rates are uniformly distributed, loss rates are 
not. However, in this application this only affects the con- 
fidence between distinguishing one photon from higher 
photon number. This occurs with very low probability, 
so this effect is not noticeable. 



V. EXAMPLE DESIGN PROCEDURE 

We now consider the details of experimental implemen- 
tation of the projective measurements from the previous 
section. We omit the A'^-port for reasons of impracticality 
[2l| and describe the optimal design of both a balanced 
TDM system and a loop architecture detector at both 
780 nm and 1550 nm wavelengths. The assumed exper- 
imental parameters for the two wavelengths are shown 
in Table [TTl using values typical of contemporary single 
mode fibre optic components. The detectors are assumed 
to be operating with 20ns gate windows and with dead 
times less than the fibre delay r. 



Component 


780 nm 


1550 nm 


Coupler loss 
Fibre loss 

Switch loss 


0.4 dB 
0.2 dB 

L = 10m (r = 50ns) 
2.0 dB 


0.5 dB 
0.8 dB 

L = 2km (r = lO^is) 
1.2 dB 


Detector 


Silicon 
Perkin-Elmer 
SPCM-AQR-13-FC 
J7dct = 60% 
Pdc = 5 X 10-'^ 


InGaAs 
id-Quantique 
id200 
Vdct = 10% 
Pdc = 9.6 X 10"* 



TABLE II: Assumed experimental parameters for the optical 
components at 780 nm and 1550 nm. 



To perform a projection onto an n photon Fock state, 
we require N > n. Thus irrespective of architecture, the 



minimum number of detection bins is N^-m = n, with 
^opt > -^min- The loop dctector exhibits the additional 
free parameter pc which may be optimised for a given 
application. Fig. [5] illustrates the optimum projection fi- 
delity achieved by both topologies at both wavelengths 
with X = 0.3. In the case of the two loop TDM lines 
we employ the coupling ratio that saturates the error 
bound from Eq. [71 where we impose a maximum error 
rate of 1%. The reason for doing this is as follows. For 
a loop TDM we always wish to minimize p^j since this 
distributes the incident field across the largest number of 
bins. However, as Pc is lowered the probability of photons 
being truncated increases. Thus, we minimize Pc subject 
to the constraint that some error bound be satisfied. In 
the case of the balanced TDM, we perform a search over 
TO, the number of stages, such that the fidelity is maxi- 
mized. The experimental parameter values found to op- 
timise the projected fidelity are listed beside each data 
point. 
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FIG. 8: Maximum achievable fidelity for the two TDM ar- 
chitectures for projection onto the Fock state shown. A^'opt = 
N-mi-a for all loop detector data points (except (1|, 780 nm 
loop). For the loop detectors the number of bins is truncated 
to A^opt . m-opt = optimum number of stages for balanced TDM 
architecture. For the loop TDM's the coupling ratios were 
chosen according to Eq. [7] so as to bound truncation error to 
1%. They were 0.99, 0.90, 0.78, 0.68 and 0.60 for = 1 ... 5 
respectively. 

Fig. [5] illustrates the corresponding pdct for each of the 
detection topologies and projections shown in Fig. [51 

A. Loop TDM detector 

Apparent from the above example is the surprising em- 
pirical result that generally A'opt = A'min. This is a con- 
sequence of dark-counts. As A^ increases the overall prob- 
ability of being affected by a dark-count increases. Thus 
it is desirable to keep A'' as small as possible, provided 
it is sufficiently large to measure n, i.e. A^min- From Fig. 
[51 it is evident that the projection fidelity decreases with 
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FIG. 9: Probability of projecting onto desired Fock state. 
|(n|i/))|^ represents the probability of finding n photons in the 
parametric downconverter output state. Note the logarithmic 
a:xes. For loop detectors the number of bins is truncated to 
No-pt- The upper line illustrates the performance of an ideal 
number resolving detector. This represents an upper bound on 
what is achievable and puts the other lines into perspective. 



photon number. The plot of in Fig. [SliUustrates 

the low probability with which these terms occur, how- 
ever. Fig. [5] also reveals a vanishingly small pdct for both 
loop detectors considered. 

Because the measurement is truncated to N time bins, 
there is some probability of photons not reaching the de- 
tector. Thus it would thus appear desirable to increase 
to improve the overall detection probability. This, 
however requires a reduction of the total loop loss (i.e. 
tftcts — > 1) to counter the effects of dark counts in order 
to maintain the projection fidelity. 

Consequently, there exist certain combinations of de- 
tector and loop parameters whereby optimum projection 
fidelities are obtained when iVopt > A'min. It is difficult to 
describe these general conditions due to their dependance 
upon both |f/)) and (n\. For example, with x = 0.3 for the 
parametric down converted state considered, N^pt = 3 for 
the 780 nm loop detector projection onto (1| (as shown in 



Fig.©. With X = 0.15 it is observed that iVopt = ^min for 
all projections. Simulated results have typically indicated 
that fidelities within 1% of the maximum attainable for 
all values of TV are achieved with iVopt = iVmin- 



B. Balanced TDM detector 

As = 2'", the minimum number of stages mmin 
required to perform a projection onto \n) is given by 
"Trnin > [1062 "1 ■ ^or thc 1550 nm balanced data plot- 
ted in Figs.[5]and[n]r7iopt — m,nin- The consequence of N 
only being able to assume powers of two is evident in the 
fidelities of the projections onto |3) and |5) appearing 
lower than the general trend would suggest. Notice that 
the projections performed by the silicon detector benefit 



from TOopt > TOmin- This is because this detector has a 
sufficiently low dark-count rate that larger values of N 
do not become corrupted by the increasing dark-count 
probability. 

Fig. [TO] indicates the combinations of pdc and 77dct 



which favour 



^opt 



for the projective measure- 



ments indicated. Note that projection onto power of two 
photon numbers places the most relaxed constraints on 
Pdc- This is a result of satisfying N^-m — n and thus the 
effects of dark counts being minimized. It is clearly seen 
that the silicon detector benefits from rriopt > rrimin in 
our application for all values of n, but the InGaAs de- 
tector does not. It should be noted that at the time of 
writing there are other InGaAs detectors available such 
as (2^ which offer better values of pdc and 77dct than our 
example, mopt > Wmin may be a worthwhile option with 
such detectors. 




10" 



40 60 



80 



100 



FIG. 10: Plot defining the detector parameters required to 
benefit from m-opt > mmin (A^opt > A^min) for projection onto 
the Fock states shown, using a balanced TDM architecture. 
X ~ 0.3 and the balanced detector component losses are as 
defined for 780 nm in Table |IT1 The characteristic detector 
parameters of pdc and r^dct are illustrated for our example 
detectors. Note the logarithmic axes. 



C. Design considerations 

Experimentalists are generally faced with the problem 
of designing a system to meet a desired specification, 
given constraints imposed by the components available. 
An example design scenario might be to determine the 
detector parameters ri^ct, Pdc required to achieve a given 
projection fidelity. Given the finite number of detector 
parameter combinations available, however, a more likely 
problem is how to maximise the fidelity for practical val- 
ues of rydot, Pdc and component losses. 

While the performance of the two architectures scale 
differently with respect to rydet and pdc, two general ob- 
servations may be made. Firstly, minimisation of pdc is 



9 



of key importance to maximise fidelity, particularly when 
projecting onto higher photon number states. Secondly, 
both schemes degrade relatively slowly with decreasing 
r^dot- Consequently in the selection of a detector, prefer- 
ence should be given to its dark noise performance. 

To select an appropriate topology, the calculations of 
Eqs. [11] and [T^] should be performed for both systems 
with component losses applicable to the wavelength of 
operation and N — N^-m initially. To determine the op- 
timum performance attainable from the loop architec- 
ture, firstly calculate the value of pc corresponding to 
the maximum desired truncation error. Secondly, iterate 
Pc through all available greater experimental coupling ra- 
tios, noting the value corresponding to greatest projec- 
tion fidelity. The above process should be repeated for 
N > iVmin when particularly low values of pdc and/or 
when significant probabilities of higher photon number 
terms exist. A choice may thus be made between the 
topologies on the grounds of projection fidelity and pdot- 

Overall, given these criteria, the best choice for a detec- 
tor to perform projective measurements is likely to be the 
balanced TDM architecture. However, both schemes offer 
useful projection fidelities and careful consideration of x, 
Vdct, Pdc and component losses is required to select and 
optimise the performance of a given architecture. Either 
architecture is suitable for the purposes of reconstruc- 
tion of photon statistics, with the loop detector having a 
potential advantage given that N may be increased ar- 
bitrarily to limits imposed by the reconstruction process 
and Pdc- In general, the presence of loss, finite detection 
efficiencies and dark counts only influence the number of 
independent measurements required to arbitrarily mini- 
mize statistical errors llal. 



photo-detection using non-number resolving detectors. 
We first considered the archetypal protocol, balanced N- 
port detection. For this architecture we considered the 
effects of loss, dark-counts and the number of bins TV. 
We also considered how to optimize this architecture un- 
der different experimental conditions. We then turned 
out attention to two forms of time-division multiplexing 
- the loop TDM, and balanced TDM. We compared these 
two schemes and how to optimize them. We also included 
results using experimentally realistic parameters. 

Our analysis has considered the major experimental 
limitations of the balanced A^-port detector cascade, and 
the loop and balanced TDM detectors. While we have 
investigated the performance of the three architectures 
from the perspective of projective measurements made 
on one mode of an optical parametric down converter, 
our treatment is sufficiently general that our conclusions 
may be extended to any projective scenario. 

The balanced TDM system appears most resilient with 
respect to non-ideal detector quantum efficiency and dark 
noise performance at the cost of requiring an additional 
non photon number resolving detector. The loop detector 
remains a useful architecture, however only for certain 
combinations of experimental parameters and projection 
operations. 
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